Quantum searches on highly symmetric graphs 
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We study scattering quantum walks on highly symmetric graphs and use the walks to solve search 
problems on these graphs. The particle making the walk resides on the edges of the graph, and at 
each time step scatters at the vertices. All of the vertices have the same scattering properties except 
for a subset of special vertices. The object of the search is to find a special vertex. A quantum 
circuit implementation of these walks is presented in which the set of special vertices is specified by 
a quantum oracle. We consider the complete graph, a complete bipartite graph, and an M-partite 
graph. In all cases, the dimension of the Hilbert space in which the time evolution of the walk takes 
place is small (between three and six), so the walks can be completely analyzed analytically. Such 
dimensional reduction is due to the fact that these graphs have large automorphism groups. We 
find the usual quadratic quantum speedups in all cases considered. 

PACS numbers: 03.67.Ac, 05.40.Fb, 42.50.Ex 



INTRODUCTION 

The theory of quantum walks describes quantum ver- 
sions of classical random walks. In these walks a quan- 
tum particle "walking" on a line, or more generally a 
graph, has different amplitudes to go in different direc- 
tions rather than different probabilities, which is the case 
for a classical walk. The time in these walks can either 
be in discrete steps [l|, Q or continuous Q . Both types 
of walks have proven to be fruitful sources of quantum 
algorithms [1, 0, H, 0] ■ A summary of both the properties 
of quantum walks and their algorithmic applications can 
be found in two recent reviews [1, Q . 

Quantum walks have been used to investigate searches 
on a number of different graphs. In these searches, one of 
the vertices is distinguished (special), and the object is to 
find that vertex. The graphs considered so far are grids 
and hypercubes of different dimensions and the complete 
graph 

In this paper, we further pursue the possibility of using 
quantum walks to speed up searches on graphs. We shall 
consider highly symmetric graphs, that is, those having 
a large automorphism group. As we shall see, on these 
graphs the quantum walk takes place in a subspace of 
small dimension, and this greatly simplifies the analy- 
sis of its behavior. The role of symmetry in quantum 
walks was studied explicitly in Ref. [l^l- The authors of 
this paper considered coined quantum walks on Cayley 
graphs of groups and found that in some cases a quan- 
tum walk on a large graph could be reduced to one on a 
much simpler graph, which they called a quotient graph. 
The situation for a search on a highly symmetric graph is 
analogous to that in the Grover search algorithm [l4| in 
which the search takes place in a two-dimensional sub- 
space. The result is that it is quite simple to find the 
behavior of the quantum search and compare it to the 



classical one. 

We shall consider searches on complete graphs and bi- 
partite graphs, and an M-partite graph, that includes 
both the complete graph and the bipartite graphs as spe- 
cial cases. We shall make use of the scattering quantum 
walk [ist . In addition, we shall show how a quantum 
search on a graph can be implemented by a quantum 
circuit containing a quantum oracle. The oracle is used 
to differentiate the properties of the special vertex from 
those of the normal vertices. 

The rest of the paper is organized as follows. We start 
in Sec. by introducing scattering quantum walks and 
show that they can be realized with a quantum circuit. 
A quantum walk search can be implemented by adding 
a quantum oracle to the circuit (Sec. ). In Sec. we 
show how symmetry of the underlying graph reduces the 
dimensionality of the problem, and then, in Sec. , we 
apply these observations to different graphs: a complete 
graph, a bipartite graph, and a special M-partite com- 
plete graph. We also show how the behavior of the search 
on the complete graph depends on the phase shift that 
the special vertex imparts to the particle upon scattering 
it. At the end of the section we also present a compari- 
son of our results with other works. We summarize our 
results in Sec. ). 



SCATTERING QUANTUM WALKS 

The scattering quantum walk, introduced in Ref. [l^, 
is a quantum walk on the edges of a graph rather than 
on its vertices. Having a graph Q — (V, E) on which the 
walk is defined, with V being the set of vertices and E 
the set of edges, the Hilbert space is defined as 



n = ^^({|m, |m, I e V, ml G E}). 



(0.1) 
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This definition states that the Hilbert space is given by 
the span of all edge states^ i.e., position states \'m^l) in- 
terpreted as a particle going from vertex m &V to vertex 
I € V, with ml £ E being an edge of graph Q. These edge 
states form an orthonormal basis of the Hilbert space, 
which we shall call the canonical basis. 

In this Hilbert space the unitary evolution is given by 
a set of local unitary evolutions defined for each vertex. 
If we set (for every m e V) = £'^{{\m, I) \l eV,mle 
E}), the set of all edge states originating on vertex to, 
and Clm = ^^({|^ m) \l GVJme E}), the set of all edge 
states ending on vertex m, then local unitary evolutions 
act as t/^"^') : Am ■ The overall unitary step opera- 

tor U acting on the system is the combined action of the 
local unitary evolutions, that is. the restriction of U to 
flm is just U^"^\ Given the initial state of the system is 
|'0init), the state after n steps is \tjjn) = U"' iV'init) and the 
probability of finding the particle (walker) in state \k,l) 
is then \ {k,l\ij)n)\'^ ■ 

Finally, let us specify our choice of local unitary step 
operators. Our vertices will be of two types, normal and 
special, and our task will be to find at least one special 
vertex. Consider a vertex I, let r(/) be the set of vertices 
connected to I by an edge, and if fc e r(Z), let V{1\ k) be 
the set of vertices connected to / by an edge but excluding 
k. The local unitary operators corresponding to both the 
normal and special vertices will act as follows: 



\Lm) 



(0.2) 



meT(l:k) 



where r^'^ and i^'-* are reflection and transmission coef- 
ficients to be chosen in such a way that is unitary. 
For normal vertices 



til) ^ 



and for special vertices, we choose 



„(0 



f« = 0. 



(0.3) 



(0.4) 



Here |r(Z)| is the degree of vertex I, i.e., the number of 
vertices in the set T{1). For both of these choices, the 
operator t/*^') is unitary, and, as we shall see, they also 
guarantee that the quantum walk has the same symmetry 
group as the graph. This choice of local unitary operators 
for the scattering walk is analogous to the choice of the 
Grover coin (see Ref. [l6j ) in a coined quantum walk. 



QUANTUM CIRCUIT SETTINGS 

Quantum search problems are often posed in terms of 
a quantum oracle. This is a black box that, when given 
an input value x, returns an output f{x). The function 
takes values and 1, and the set of inputs for which it 
takes the value 1 is typically small. The objective is to 



IV'init) 



|0) 



Wo 



|0) 



FIG. 1: A logical circuit (network) that implements a single 
step of a scattering quantum walk search, which makes use 
of the quantum oracle CVf. The first input corresponds to a 
quantum walker originally prepared in the state It/"!!!!!). The 
second input represents a vertex state, while the third input 
represent an ancillary qubit. 



find at least one input for which the function is equal to 
1. What we would like to do in this section, is provide a 
connection between this phrasing of a search problem and 
a search whose object is to find a set of special vertices 
on a graph. We shall accomplish this by showing how an 
oracle can be used in a quantum circuit to implement a 
quantum search on a graph. 
Our oracle has the action 



CVf \k) ® |to) = |fc) ® \m®J{k)) 



(0.5) 



where k is the input and / is the function that determines 
the action of the oracle. We shall assume that /(fc) is 
either or 1, though this restriction is not necessary. The 
second system (ancilla) is, therefore, a qubit, and © is 
addition modulo 2. In addition to this operator, we need 
one that will link the values of the function to the action 
of local unitary evolutions at the vertices. In particular, 
if Hs a vertex of the graph, the f{l) = will correspond 
to a normal vertex and /(/) = 1 will correspond to a 
special vertex, i.e., the type of vertex we are trying to 
find. If Z is a normal vertex, the local unitary operator 
corresponding to it will be denoted by Uq \ and if it is a 
special vertex the local unitary will be denoted by Ui\ 
Our quantum circuit will act on a tensor product of the 
Hilbert space Ti, for the quantum walk [sec Eq. (|0.ip ] , the 
Hilbert space for vertices, 7i„, and a qubit Hilbert space 
7^2- The vertex space is given by 



n.^e{{\i)\i&v}). 



(0.6) 



We now define an operator, CUf, acting on this space in 
the following way 

CUf{\k, I) ® \l) ® \c)) = (u[^^ \k, O) ® \l) ® \c) . (0.7) 

This equation does not completely specify the actions of 
CUf. In particular, it does not specify its action on states 
of the form \k, I) (E) \l') (g) |c), where I ^ I', but we will only 
need to consider its action on states of the form given in 
the previous equation. 
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The quantum circuit that implements one step of our 
quantum walk search is given in Fig. [TJ The first in- 
put stands for a state of quantum walker (i.e. any su- 
perposition of edge states), the second for a vertex state, 
and the third for an ancillary qubit. The input state is 
iV'o) = iV'init) (81 |0) (g) |0) where IV'init) is a general state 
in the edge Hilbert space 7i, 



|'0init) = ^ a^i\m,l), ^ \a„ii\ 



1. 



The state |0) in the second input in Fig. [T] is one of the 
vertex states, which, besides labeling a particular vertex, 
will also serve as a reference state. First, we apply the 
operator CW2 which maps the state |m, I) ® |0) in H<^Hv 
to \m,l)®\l). Such a unitary can be implemented, e.g. as 
presented in Refs. 17, 3| and references therein. After 
this operator is applied, the initial state becomes 

i^jo) ^ 1"^' ^) ® 1^) ^ 1^^ = 1^1^ • 

Next, we apply the quantum oracle given by Eq. (jO.Sp to 
the vertex state and the qubit, yielding 

1^1) a,ni \m, I) (g) 10 (g) 1/(0) = 1^-2) • 

mlGE 

Now we can apply the CUf operator from Eq. (|0.7p to 
the state, which gives us 

1^2) ^ Y (^/S) l'"'^)) ® I') ® = 1^3) ■ 

ml£E 

Here it is clear that both ancillary systems act as controls, 
and they will have to be reset before we can use the circuit 
again. That is the task of the remaining two gates in the 
circuit. Now, because the local unitary operator ^^j';) 
acts only on the edge space, the state can be rewritten 
as 



1^-3 



InieE 



|/,m)®|0® 1/(0)- 



Application of the quantum oracle again resets the qubit 
state to |0). However, to reset the vertex state, we can- 
not use CW2 operation as before, since the information 
contained in the edge state about the vertex state has 
moved from the second to the first slot. Therefore, we 
define the opcratorCVKi, which maps \l,m)^\l) m'H<E)'Hv 
to \l,m) (g) |0), and apply it to our state, giving 



Y bi,n\l,m)(E)\0)(E)\0). 



lm£E 



We have thus performed one step of the walk and reset 
the ancillas, so that the circuit can be applied again to 
perform additional steps of the walk. 



SYMMETRY CONSIDERATIONS 

For simplicity, suppose we have a graph Q = (V, E) 
with only one special vertex, and let A be the group of 
automorphisms of the graph that leave the special vertex 
fixed. An automorphism a of ^ is a mapping a :V ^ V 
such that for any vertices wi , W2 £ V there is an edge 
connecting a{vi) and a{v2) if and only if there is an edge 
connecting vi and W2- Each automorphism, a, induces a 
unitary mapping Ua on the Hilbert space of the graph Q 
given by Eq. (jO.ip . such that Ua\vi,V2) = |a(wi), 0(112)). 
Suppose now that Ti, can be decomposed into m sub- 
spaces. 



where each Tij is the span of some subset Bj of the canon- 
ical basis elements and is invariant under Ua for all a ^ A. 
We shall also assume that each Tij does not contain any 
smaller invariant subspaces m. 

Next, in each invariant subspace we can form a vector 
that is the sum of all of the canonical basis elements in 
the subspace. 



1 



E 



\Vl,V2) 



(0.8) 



where dj is the dimension of Jij. This vector satisfies 



UAwj 



for all a <E A. Moreover, it is the only 



vector in Tij that satisfies this condition. Define S ~ 
= l,2,...,m}), and note that S = {\'ip) G 
T-L\tJa\i^) = Va e A\, and that the dimension of S is 
simply the number of invariant subspaces. Now suppose 
that [?7, C/q] = for all a A. This implies that if 
ij^\'4,) = then Ua(j\^^) = and if IV') G 5, then 

G S. Therefore, the subspace S is closed under the 
action of the step operator U . Correspondingly, if the 
initial state of the walk is in 5, then we only need to 
consider states in S to describe the state of the walk at 
any time. If the automorphism group is large, then S 
can have a much smaller dimension than Ti. 

Now let us demonstrate that the unitary operator J7, 
defined by the local unitary operators in Eq. (|0.2p . does, 
in fact, commute with all of the automorphisms of a 
graph that leave the special vertex fixed. If these op- 
erators commute when applied to all of the elements of 
the canonical basis, then they commute. As before, r(w) 
is the set of vertices in V that are connected to the ver- 
tex w, and, if v' G r(w) then r{v;v') = T{v) — {v'} and, 
finally, |r(w)| is the number of elements in r(7;). Then 
we have that 



UaU\vi,V2) = -r(''=)|a(w2),a(«i)) 



E 
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We also have 



UUa\vi,V2) 



v^T{a(v2);a{vi)) 



First note that the reflection and transmission ampli- 
tudes in this equation are the same as those in the previ- 
ous equation, i.e., r^''^) = A^M] and ^^^^^ = ^I''^^^)] — 
this is a consequence of the fact that |r(z;2;t'i)| = 
|r(a(w2); a(wi))| and the fact that the special vertex is 
mapped into itself by the automorphism. We also have 
that T{a{v2)',a(yi)) = {a(w)| v £ T(v2]vi)\, so that 
the sums in the two equations are identical. Therefore, 

Ua(j\vi,V2) = tjUa\vi,V2). 

Somewhat similar considerations appear in Ref. [31, 
but the situation there is complicated by the fact that the 
analysis was performed for a coined quantum walk. This 
requires the edges to be colored, and the automorphisms 
to respect the coloring, that is, edges must be mapped to 
edges of the same color. The situation for the scattering 
quantum walk is much simpler, because no edge coloring 
is required. 



EXAMPLES OF SEARCHES 
Search on a complete graph 

Let us consider a complete graph with TV vertices 
(see Fig. [2]) . Each vertex of this graph is connected to 
all of the other vertices by an edge, so the graph has 
N{N — l)/2 edges which define the Hilbert space of di- 
mension dimTY = N{N — 1). Let v be the number of 
special vertices. Then, without a loss of generality, we 
may label them as k = l,2,...,z; and corresponding lo- 
cal unitary evolution is defined by Eqs. (|0.2p and (|0.4p . 
Normal vertices are labeled as k = u + 1, w + 2, . . . , A^, 
and have local unitary evolution defined by Eqs. ()0.2p 
and ()0.3p . Transmission and reflection coefficients for all 
normal vertices a are the same, with |r(a)| = — 1. 

The Hilbert space for the walk defined on this graph 
can be decomposed into four subspaces characterized by 
four vectors according to Eq. (|0.8p . 



1 



N 



\Wi) 



\W2) 



V N 



1 



Jv(N - v) 1 , , 1 

V \ / a=l b=v+l 



E E 



1 



N N 
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ViN-v)iN-v^l) ^^^^ 




FIG. 2: An example of a complete graph with N — 7 vertices 
out of which V — 2 are special (white ones). A solution for 
a scattering-quantum-walk search on such a graph leads to a 
reduction in dimensionality of the problem to four dimensions. 



W4 



^ V V 

=T^TEEK^) 

v{v - 1 r— ^ 



a=l b=l 
b^a 



These are equal superpositions of all edge states directed 
from normal vertices to special vertices, from special ver- 
tices to normal vertices, connecting normal vertices and 
connecting special vertices, respectively. Unitary evolu- 
tion, as was shown in Sec. , can be described by the 
evolution within the subspace S spanned by four vectors 
\wk), k = 1,2, 3, 4, and is described by a 4 x 4 matrix 



U = 



/ q s 0\ 

e'"* 

s -q 

V c'-^y 



(0.9) 



where 



-r + t{v - 1) 



2v 



s = yjl - (j2 = t^Jv{N -u-l). 

Note that the subspace spanned by the vectors \'Wk),k = 
1,2,3, is decoupled from the subspace spanned by the 
vector |w4). For the initial state we choose an equal su- 
perposition of all edge states, i.e.. 



IV^init) = 



1 



^N{N-l) 



E i^'O 



k,l)£H 



This initial state is a natural choice, since no set of edge 
states is preferred. We have that IV'init) € S and can 
be written as a superposition of the states \wk) for k = 
1,2,3,4, 



init/ 



N{N ~ 1) 



{\WI) + \W2)) 



'{N ~v){N -v-l) 
N{N - 1) 




N{N - 1) 
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We shall analyze what happens for two choices of a 
phase shift, namely, = and cj) = tt. For <j> = 0, we find 
that the initial state can be written as a superposition 
of eigenstates of U with eigenvalue equal to unity. In 
particular, the following two states are eigenstates of tl 
with the eigenvalue equal to unity: 

/ 1 \ 
1 



"0 



1 



3 + q 



1-9 



1 + 9 

V J 

14) = (0,0,0,1)^, 

and the initial state can be expressed in terms of them 
as 



IV'init) 



'{N~v){N + v-l). . 



N{N - 1) 



NiN-iy 

This, however, means that the walk goes nowhere; the 
particle remains in the initial state for the entire walk. 
In this case, the quantum walk gives us no advantage 
over a classical search. Classical searches on edges are 
discussed in Appendix A. 

The value of = tt is much more interesting. The state 
of the walk after n steps, derived by using the formula 



IV'n)=$^(A"(uA|7^init»kA> 



(0.10) 



where \ux) is the eigenvector corresponding to eigenvalue 
A of unitary operator U, is, in the limit N ^ v > 1, 



init/ 



^y2v{N - l)sin(2n + 1); 



V2w(^- l)sin(2n - 1)- 







1) cos n9 



where 



tan 6 = 



N -v-1 ' 



(0.11) 



(0.12) 



We see that the probability amplitudes for edge states not 
connected to special vertices are approximately equal to 
zero when On = n /2. If we measure the position of the 
particle after this many steps, it will, with probability 
close to unity, be located on an edge connected to one 
of the special vertices. This implies that the number of 
steps needed to find one of the special vertices is of the 
order 0{\jN/v) for large A'' which is a quadratic speedup 
over the classical algorithm that needs 0{N/v) steps to 
do the task when searching an unstructured database. 



It should be noted that for the case v = 1 there is no 
vector \wi) and the dimension of the problem is reduced 
to three, yet the results remain. 



Role of the phase shift 

As we have seen in the previous section, the differ- 
ences in the behaviors of quantum walks with a phase 
shift and a phase shift tt are substantial. To see how 
this change occurs we will present numerical results for 
< (/) < 27r for the case of a single special vertex. The 
behavior of the walk is depicted in Fig. [3l where we see, 
that the further we arc from = tt, the smaller is the 
probability P of finding the particle on an edge that is 
connected to the special vertex. 

We can see from Fig. |3] that P is symmetric about 
(j) = TT. This property can be shown in the following way. 
If we make the substitution — > 27r — (/) in the unitary 
U from Eq. we find that (j changes to U*. This 

result together with the fact, that the initial state has 
real-valued coefficients when expanded in the canonical 
basis, yields 



U" 



So the resulting components arc now the complex con- 
jugates of the original; however, the probabilities remain 
the same. 

The value = tt is special, because for any other 0, 
the probability P never reaches unity. The significance 
of = TT leads to the following problem: Suppose we are 
given a complete graph having one special vertex with 
an arbitrary but known phase shift. Our task is to find 
the special vertex. In order to find it efficiently, we need 
to wait for the optimal number of steps before perform- 
ing the measurement. If we were given only one chance 
to measure, we would have to wait until the probability 
of finding a particle on one of the edges connected to 
the special vertex reaches its maximum. However, if we 
arc able to repeat the experiment an arbitrary number 
of times, we may, after an unsuccessful search, do the 
experiment again. In that case, the optimal number of 
steps before measuring is different — it is the number of 
steps, for which the average number of steps is minimal. 

Let (m) be the probability of finding the particle on 
an edge connected to the special vertex after one repeti- 
tion of the experiment assuming that to steps of the walk 
are taken before the measurement, and the phase shift is 
(j). Then the average number of steps n^,™ to be taken 
when measuring after m steps on a graph with phase shift 
(j) is given by 



= y^^P^,m{k)kr 



k=l 
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FIG. 3: (Color online) The density plot of the probability P, 
taken for a quantum walk on complete graph with A'" = 256 
and V — 1, shows the convergence of the side ridges to the 
maximum (dark gray) for (j) = -k with the probability P = 1 
(blue dots). White areas correspond to the minimal probabil- 
ity of finding the particle on an edge connected to the special 
vertex. The thick red line represents the optimal number of 
steps that are to be taken before the measurement so that the 
search algorithm for the special vertex would require minimal 
average number of steps. 



where 



P4>..mik) - [1 - P0(m)]^-ip^(m) 



is the probability of finding the particle on an edge con- 
nected to the special vertex after k repetitions of the 
experiment (and no sooner). After a short evaluation we 
find 
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FIG. 4: Average number of walking steps as a function of 
the phase shift (j). This plot illustrates the efficiency of 
algorithms — the classical blind search (dashed), the classical 
search with memory (dotted), and the quantum search (thick 
solid). We consider a graph with N — 256 vertices. We 
see that the minimal value of the average number of steps is 
achieved for the quantum search with = tt. 



Bipartite graph 

Let us now consider a bipartite graph consisting of 
two sets of vertices in which each vertex in each set is 
connected to every vertex in the other set. but there is 
no connection between vertices in the same set. Let us 
further suppose there are Ni vertices in the set 1 and N2 
in the set 2. Of all vertices in the set 1, Vi are special 
vertices and pi = A^i — vi are normal vertices, and in the 
set 2, V2 are special vertices and p2 = N2 — V2 are normal 
vertices — see Fig. [S] The action of the unitary operator, 
U acting on a state entering a normal vertex in the set 1 
is given by Eq. ([(T^ with(2a| 



2 

iV2^ 



ri = 



N2-2 
No 



m 
P0(m) 

The optimal number of steps, mopt, in each experiment 
is given by the value of m that minimizes n^,m- We will 
denote the average number of steps to be taken in the 
optimal case as n(<j)) = n^^„iapt ■ This is plotted in Fig. U) 
For = TT, the quantum approach is faster than the 
classical; it has a quadratic speedup, even in the nonop- 
timal search via the maximum probability method (see 
Fig. HI). However, for the case of ^ we see that the op- 
timal number of steps before measuring is — we do not 
evolve the system, we just measure the initial state — and 
the average number of steps to find the special vertex 
reaches a value close to N. This is the same situation 
we would have if we performed a classical blind search 
in which wc randomly pick vertices at each step without 
remembering the past choices. 



If it acts on a state entering a normal vertex in the set 2 
its action is again given by Eq. (|0.2p but with 

2 TVi - 2 

The action of U on a state entering a special vertex is 
given by Eqs. (|0.2[) and (|0.4p . We shall number the ver- 
tices in the set 1 as 1 through A^i and those in the set 
2 as iVi + 1 through A^i + N2. In analyzing the walk, 
we can, without loss of generality, assume that vertices 1 
through vi in the set 1 are special, and vertices A^i + 1 
through A^i -|- V2 are special as depicted on Fig. [S] 

There are now eight invariant subspaces, two of which, 
consisting of edges that connect the sets of special ver- 
tices, are decoupled from the rest. So our problem is es- 
sentially six-dimensional. We can see this by first defining 



7 




«1 +Pl+ t-'2 + P2 



«1 +Pl 



FIG. 5: General bipartite graph having vi special (black cir- 
cles) and pi normal (white circles) vertices in the first set 
and V2 special and p2 normal vertices in the second set. This 
problem can be reduced to an eight-dimensional one with a 
two-dimensional decoupled subspace. The remaining six di- 
mensions can be reduced to three by performing two steps at 
a time due to the oscillatory behavior of the walk in bipartite 
graphs. 



the vectors 



|Woi) = 



\W02) = 



N1+V2 
E 



\k,j) 



j=l k=Ni + l 
vi N1+V2 

E E b-'^)- 

j=l k=Ni + l 



The vector \woi) is a superposition of states leaving spe- 
cial vertices in the set 2 and entering special vertices in 
the set 1, and |wo2) is a superposition of edge states leav- 
ing special vertices in the set 1 and entering special ver- 
tices in the set 2. We have that U\wni) = exTp{i(j))\wo2) 
and U\wq2) = cicp{i(j))\woi) ■, so that these vectors decou- 
ple from the rest of the problem. 
Now let us define 



1 



\Wii) = 



\Wl2) 



VI N1+N2 

, y y \k,j) 



JVi+tl2 



1^13/ 



y/PlP2 



E E 

Ni 

E 



N1+N2 

E 



\k,j) 



These vectors consist of different sets of edge states en- 
tering the set 1. Let us denote their span by Si. Let us 
define the vectors 



1 



JV1+7V2 



\W21/ 



\W22) 



\W23) 



E E b'^)' 



j = l k=Ni+V2 + l 



Ni 

E 



N1+V2 



Ni N1+N2 

y y 

^ 3=vi + l k=Ni+V2 + l 



.hk). 



These vectors consist of different sets of edge states en- 
tering the set 2. We shall denote their span by 1S2. 

Finally, let us consider the action of the unitary oper- 
ator that advances the walk one step on these states. We 
find that 



U\wn) 

U\W12) 



6*1^21), 
[tiiv2 - 1) 
[h{P2-l) 



and 



U\w2l) 
U\W22) 
U\W23) 



MVI - 1) 
[t2{p2-l) 



ri]\w22) + tiy/V2P2\w23), 
ri]\w23) + tly/V2P2\w22), 

f2]|wil) + t2^/vrPl\Wl3) , 
f2]|wi3> + t2^/ViPi\wii). 



From these equations, wc see that U maps Si into 52 and 
vice versa. This also implies that C/^ maps Si into itself 
and 1S2 into itself. Therefore, if we consider a walk with 
an even number of steps, our six-dimensional problem 
turns into two three-dimensional ones. This is what we 
shall do. 

Wc shall look in detail at what happens in ^2 ; the case 
of Si is similar. We shall only consider the case = tt. 
First, let us define the quantities 

qi = -ri + ti{v2 - 1), Si = tiy/V2P2, 
q2 = -r2 + t2{vi - 1), S2 = t2y/ViPi, 

and note that -I- s| = 1 for j = 1, 2. If we denote 
restricted to ^2 by M, then the matrix for M in the basis 

\w2i), \w22), and |w23) reads 

-q2 -S2 

M = I S1S2 -qi -<?2Si 

-qiS2 -si qiq2 

In finding the eigenvalues and eigenvectors of this matrix, 
we assume that the number of special vertices is small, 
that is, (vi/Ni) <^ 1 and {V2/N2) <?; 1, and we further 
neglect higher-order terms. In addition, we set xi = 
2vi/Ni and X2 = 2V2/N2. 

For the initial state of our walk, wc will choose the state 
that is an equal superposition of all edge states entering 
the set 2. The behavior of the walk is similar for any 
initial state that is a linear combination of this state and 
the state that is an equal superposition of all edge states 
entering the set 1. Note that the state that is an equal 
superposition of all edge states in the graph is of this 
form. For simplicity we choose to start the walk in the 
equal superposition of all edge states entering the set 2. 
In this case the entire walk, to very good approximation, 
takes place in the subspace 52. Our initial state can be 
expressed as 



IV'init) 



1 



+ \/V2Pl\w22 



{y/VlV2\wo2) + y^VlP2\w2l) 
VP1P2IW23)). 
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This vector is not entirely in 1S2 due to the |w;o2) compo- 
nent. But the component that is not in 1S2 is small, and 
it stays small throughout the evolution. Neglecting this 
small component, wc find that 

f- y'xilixi + X2 ) sin(2n6l)\ 
[/^"IV'init) ^ Vx2/(.Ti+X2)sin(2n0) . 

V - cos(2n6') / 

From this equation, we see that, when 

7r 7r 

2"=7T7T = — , (0-13) 
20 2^2(^1 +X2)' ^ ^ 

we are with certainty on an edge connected to a special 
vertex. If after this many steps we measure the particle to 
determine which edge it is on, with probability xi/{xi + 
X2) we find it on an edge connected to a special vertex 
in the set 1, and with probability x^jixx + X2) we find it 
on an edge connected to a special vertex in the set 2. 

In order to get a better understanding for this solution, 
let us consider the case v\ = V2 = 1. In this case 

When 2n = the probability of finding the parti- 

cle on an edge connected to the special vertex in the 
set 1 is N2/{Ni + N2) and the probability of being on 
an edge connected to the special vertex in the set 2 is 
N'i/{Ni + N2). Now let us suppose that N2 > Ni and 
ask the question: How many steps would it take to find 
the special vertices in each of the sets? We find that the 
number of steps in the walk is tt/29 ^ V-^i- The num- 
ber of times the walk would have to be repeated in order 
to find the special vertex in the set 1 is 0(1), while the 
number of times to find the special vertex in the set 2 is 
0{N2/Ni). Therefore, the total number of steps to find 
the special vertex in the set 1 is O(V-^i), and the total 
number of steps to find the special vertex in the set 2 
is 0{N2/\^7^)- A classical search would require 0{Ni) 
steps to find the special vertex in the set 1 and 0(A^2) 
steps to find the special vertex in the set 2. 

Another interesting case is when all of the special ver- 
tices are only in one set. For example, if they arc all in 
the set 1, then X2 = 0, and we see from Eq. (j0.13p that 
the number of steps that concentrates the probability on 
edges connected to special vertices is 0{^J Ni/vi), which 
does not depend on the number of vertices in the set 2. 
Therefore, if the special vertices are only in one set, the 
size of the other set does not affect the number of steps 
necessary to find a special vertex. 




FIG. 6: A/-partite complete graph consists of M sets of ver- 
tices where each set contains vertices. For every pair of 
vertices not belonging to the same set there exists an edge, 
while there is no edge connecting any two vertices within the 
same set. In one of the sets we consider one vertex to be 
special. 

vertices. Vertices within sets are not connected while 
for any pair of vertices not belonging to the same set 
there exists an edge connecting them (see Fig. E]). We 
see that by taking M = 2 we obtain the bipartite graph 
discussed in Sec. . However, if we suppose that there is 
only one vertex in each set, we get a complete graph, 
which was discussed in Sec. . In both of these cases, a 
quadratic speedup of classical searches has been demon- 
strated. For the graph we are presently considering, M is 
its chromatic number. The chromatic number of a graph 
is the number of colors one would need to color its ver- 
tices so that no two vertices connected by an edge have 
the same color. 

To simplify the analysis, let us suppose, that there is 
only one special vertex and both M and N are large 
(see Fig. |6]). We will label the vertices by {m,n), where 
m = 1,2, ... ,M is the number labehng a group (set), 
and n is the number labeling vertices within the group. 
We shall further assume that the special vertex is (1,1). 
It turns out that this problem is five-dimensional, where 
the subspace S in which the walk is happening is defined 
by five vectors according to Eq. (|0.8p : \wi) consisting 
of edges entering the special vertex, \w2) consisting of 
edges leaving the special vertex, jwa) and |u;4) consisting 
of the edges leaving and entering normal vertices in the 
set 1, respectively, and {w^) consisting of all of the edges 
connecting vertices in sets 2 through M. Note that in 
the case M = 2 we do not have the vector jws). 

The unitary operator U advancing the walk by one step 
is again defined by Eqs. (|0.2|) and (|0.4|1 for the special 
vertex and by Eqs. ()0.2p and ()0.3|) for normal vertices 
with 



M-partite complete graph 



As a generalization of previous cases we may consider 
an M-partite complete graph consisting of M sets of N 



By analyzing the action of U on vectors \wk) for 
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fc = 1, 2, . . . , 5 (see Appendix B), neglecting higher-order 
terms, and, finaUy, employing Eq. (|0.10|) we find that 



IVinit) 



/ sin (j>n \ 
— sin (f>n 



\\/2 cos (j)n / 



where cos(f> ~ 1 — 1/MN. We observe, that if <j>n = ^, 
only the first two components are nonzero, leading to the 
probability of finding the sought element being close to 1. 
In this case n ~ \/MN , which is again quadratic speedup 
over classical search, in which 0{MN) steps would be 
needed. 



Comparison with other works 

There are several papers (to the authors' knowledge) 
that have also studied searches on graphs. We would like 
to compare the results in these papers with our results. 

The problem discussed here is, in fact, similar to the 
one in Ref. 



I4I , where a search on a hypercube was stud- 



ied. There an algorithm for finding the (only one) special 
vertex with a probability of approximately 1/2 was pre- 
sented, which is in contrast with our finding of a probabil- 
ity close to unity. In our case the probability is (almost) 
equally split between two possible sets of edge states, 
those leaving the special vertex and those entering the 
special vertex. In the case of Ref. [l2| , this corresponds 
to either finding the particle on the special vertex (with 
an arbitrary coin state) or finding it on one of the neigh- 
boring vertices with the coin pointing to the special one. 
The authors of Ref. discard the latter type of states 
as an unsuccessful search. Similar findings are given in 
Ref. [111. 

In Ref. [lo| a quantum walk that performs an exact 
Grover search is discussed. The graph the quantum walk 
is performed on is a complete graph that has loops added 
to each vertex. Coins (equivalents of our local unitary 
evolutions) for normal vertices are Grover coins. How- 
ever, coins for the special vertices are "minus" Grover 
coins. This leads to almost the same evolution as in the 
case of a Grover search; the only difference is that one 
step of the Grover algorithm corresponds to two steps of 
the quantum walk. 

Finally, let us summarize the situation as it presently 
stands. Quantum- walk searches have two figures of 
merit, the number of steps necessary to find a special ver- 
tex and the probability of finding it after a specific num- 
ber of steps have been made. On a complete graph with- 
out loops [the result obtained in this paper in Eq. (jO.lip ] 
one needs v^-times as many steps as in the Grover search 
for the corresponding problem (a search on N objects), 
and after this many number of steps the probability of 



finding the special vertex is equal to unity. On a com- 
plete graph with loops twice as many steps as in 
the corresponding Grover search are required, and the 
probability of finding the special vertex is again equal to 
unity. A rigorous comparison of these properties on a 
hypercube can be found in Ref. [ll| , where adding loops 
to the graph incurs the cost of increasing the number of 
steps in the walk. 



CONCLUSION 

We have considered scattering-quantum- walk searches 
on several examples of highly symmetric (complete, bi- 
partite, and M-partite) graphs where some of the vertices 
are special. We have shown how an oracle in a quantum 
circuit can be used to mark these special vertices. The 
symmetry of these graphs leads to a significant reduction 
in the dimensionality of the problem. For all of the types 
of graphs we considered, we found a quadratic speedup 
over a classical search. These results were obtained by 
taking the phase shift of special vertices to be tt. Taking 
0=0 results in a trivial evolution, with constant prob- 
abilities of finding the particle in any edge state. In this 
case quantum walks reduce to the classical blind search 
on edges. That is, the quantum search is reduced to its 
classical counterpart. 

We have also studied the change in behavior when 
changing the phase shift (p for a complete graph with one 
special vertex. While for (j) — tt wc have a nontrivial be- 
havior suitable for searches in these graphs, for </) = we 
get only a static walk. Cases between these values were 
explored numerically. As a measure for suitable compar- 
ison between different choices of 0, we chose the average 
number of steps that need to be taken to successfully find 
the special vertex. We note that the quantum algorithm 
is always at least as fast as its classical counterpart. 

Finally, we have made comparisons of our results with 
other works, putting this work into perspective. It is 
hoped that the approach used here makes the evaluation 
of these cases easier and more simple to interpret. 
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CLASSICAL SEARCHES 

When employing quantum oracles in searches for an 
element in an unstructured database, the classical search 
may be viewed as a quantum query not capable of using 
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superpositions. In this way one may consider two types 
of classical searches — blind searches and searches with a 
memory. We present both types of search below. 

Blind searches 

Classical blind searches are searches where previously 
chosen elements are not marked, so that they may be 
picked up again. This means that the probability P of 
choosing a special clement remains the same after every 
unsuccessful step of the search. The probability of finding 
a special element after k steps, hence, is 

Pfe = (1 - Pj^'-ip, 

which is the probability of not finding the element in A: — 1 
steps and finding it at the k-th step. The average number 
of steps we need to take in order to find a marked element 
is expressed as usual: 

oo 

n = ^Pfcfc. 

fe=i 

In the case when the oracle function marks v vertices out 
of iV, we have P = v/N so that 

1 N 
P V 

Classical searches with memory 

Let us suppose that we mark previously chosen el- 
ements, and do not choose them again in subsequent 
steps during the search. The probabilities Pk of find- 
ing a marked element now change for each new element 
chosen: 

-V N -v~l 



Pk 



N 

N - 



N - 
- k- 



N -k 
{N-v)\v 
iV! 



{N 



N -k- 

-ky. 



{N - V ~ k + 1)1' 

This is a product of (fc — 1) probabilities of not finding 
a special vertex in successive searches, where after each 
unsuccessful search we remove the selected normal vertex 
from the search. In this case, the average number of steps 
taken in order to find a marked vertex is 

N-v+l 
k=l 

After some evaluation this yields 

iV + 1 N 
n = ~ — . 

V + 1 V 

In both cases we may observe that the efficiency remains 
the same, of the order N/v. 



M-PARTITE GRAPH 

The explicit expressions for the vectors . . . , \wz) 
read 



M N 



\W2) 



[N{M - 1)] 



\ ^)11/2 5Z £ |K7l),(l,l)), 



m— 2 n— 1 

M N 



[N{M - 1)] 



^TWI £ £ |(l,l),(m,n)). 



in—2 n—1 



1 



[iV(Ar-l)(A/-l)]i/2 

M N N 

^ 5Z H |(l,ni),(m,n2)), 

m— 2 ni— 2 n2 — 1 

1 



W4 



[iV(Ar- 1)(M- l)]i/2 

M N N 

^ H 51 H |(m,n2),(l,ni)), 



m— 2 ni— 2 71-2 — 1 

1 



M M 



W5 



[iV(M- l)(A/-2)]i/ 



1/2 X! X! 



mi— 2 m2— 2 



N N 



m — 1 712—1 

The actions of U on these vectors are 
U\w{) = e'^\w2), 

U\w2) = -r\wi) + t^W^\wi) + t^N{M - 2)1^5), 

+t[N{N^l){M-2)Y'^\w,), 
U\wi) = |u)3), 



U\wz) = t^N{M~2)\wi) 

M — 

+t[N{N - l)[M - 2)Y/^\w,) + ^fj^ks). 

These equations give us a 5 x 5 matrix, whose eigenvalues 
in the large-M and large- limit are approximately — 1, 
±i — (1/Af), and 1 ± zi/t, respectively. Only the eigen- 
vectors corresponding to the eigenvalues 1 ± iVf have a 
significant overlap with the initial state, so it is only these 
eigenvalues that determine the behavior of the walk. 
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